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In this paper we study the isotropic cases of static charged fluid spheres in general rela- 
tivity. For this purpose we consider two different specialization and under these we solve 
the Einstein-Maxwell field equations in isotropic coordinates. The analytical solutions 
thus we obtained are matched to the exterior Reissner-Nordstrom solutions which con- 
cern with the values for the metric coefficients and e*^ . We derive the pressure, density, 
pressure-to-density ratio at the centre of the charged fluid sphere and boundary R of the 
star. Our conclusion is that static charged fluid spheres provide a good connection to 
compact stars. 
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1. Introduction 

The topic of static charged fluid spheres in general relativity is a challenging 
issue and has given rise to interesting studies, even if their number is not so high. 
Considering the research of the interior of stars in connection to their late stage 
evolution when the general relativistic effects play an important role, we notice the 
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interesting work of TolmanlHthat yields a class of solutions for the static, spherically 
symmetric equilibrium fluid distribution. This important result was followed by 
some generalizations made by Wyman ^ LeibovitzEl and Whitman SI. After them, 
Bayin used the idea of the method of quadratures and gave new astrophysical 
solutions for the static fluid spher es. 

In recent years, Ray and Das ISEE! ^nd Ray et. al. — have developed some in- 
teresting solutions for the static charged fluid spheres in general relativity following 
the above line of thinking. Ray and Das ^ performed the charged generalization 
of Bayin's work ^ related to Tolman's type astrophysically interesting aspects of 
stellar structure ^. In this light, they actually considered that even in the case of 
stellar astrophysics there are physical irnplications of the Einstcin-Maxwcll field 
equations. In other papers Ray and Das performed a study of some previous 
solutions considering the phenomenological connection of the gravitational field to 
the electromagnetic field, and demonstrated the purely electromagnetic origin of the 
charged relativistic stars given by TolmanlHand Bayin IS. The existence of this type 
of a stro physical solutions is thought to be a probable extension of Lorentz's conjec- 
ture^I^that electron-like extended charged particle possesses only 'electromagnetic 
mass' and no 'material mass'. In their latest work Ray et al. ^ have considered 
Tolman-Bayin type static charged fluid spheres in general relativity and they have 
found out many interesting results, specially the cases which give support to the 
charged spherical models in connection to normal stars. It is argued by them 1^1 that 
due to the inclusion of charge and by a suitable choice of charge part, pressure and 
density function could be a decreasing function of radius from centre to surface in 
contrary to Bayin's case^. 

In this connection we also notice other two works, the first one elaborated by 
Varela IHJ, which presents a neutral perfect fluid core bounded by a charged thin 
shell, and which demonstrates that it is possible to construct extended Reissner- 
Nordstrom sources with everywhere positive mass density, classical electron radius, 
electromagnetic mass, an d ev erywhere non-negative gravitational mass. The other 
work was done by Ivanov that studied the interior perfect fluid solutions for the 
Reissner-Nordstrom metric using a new classification scheme. In addition, he found 
general formulae in some particular cases, presented explicit new global solutions 
and made a revision of some kno wn solutions. In connection to the singularity 
problem it is argued by Ivanov that the presence of the charge function serves 
as a safety valve, which absorbs much of the fine-tun ing nec essary in the uncharged 
case. However, it is also believed by some authors l ^'^ l ^'^ l that in the presence of 
charge the gravitational collapse of a spherically symmetric distribution of matter 
to a point singularity may be avoided through counterbalancing of the gravitational 
attraction by the repulsive Coulombian force in addition to the thermal pressure 
gradient due to fluid. 

In continuation of the above theme, especially that of Ray et al. in the present 
article we have tried to solve the Einstein-Maxwell field equations in isotropic coor- 
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dinate system and derived expressions for pressure and density. We have found out 
conditions for the boundary of the charged sphere. The exterior Reissner-Nordstrom 
soliition is compared and constants of integrations are expressed in terms of mass 
and radius. The mass-radius and mass-charge relations have been found out for 
various cases of the charged matter distribution along with the pressure-to-density 
ratio at the centre of the charged sphere. 

Our paper is organized as follows: in Section 2 we introduce the Einstein-Maxwell 
field equations for the static charged fluid spheres in general relativity, we determine 
their analytical solutions making some assumption for the metric coefScients e'' and 
e** and we consider two particular cases. In Section 3 we find out boundary condi- 
tions, which consist in a vanishing value for the pressure and a specific value for the 
boundary R of the star. The expression for the radius of the star R is established 
in Section 4 that is dedicated to a detailed analysis of the obtained solutions. The 
role of some parameters and integration constants in the evolution of pressure p, 
density p and radius of the star R is discussed, with emphasis on some particular 
values. In Conclusions we enlighten the physical significance of the results. 



2. The Einstein-Meixwell field equations and their analytical 
solutions 

The static spherically symmetric matter distribution corresponding to the 
isotropic line element is given by 

= e'di^ - ei'idr^ + r'^dO^ + r'^sin'^edci)''), (1) 

where y and are functions of the radial coordinate r only. 

The Einstein-Maxwell field equations for the above line element can be written 

as 

where the total charge with a sphere of radius r in terms of the 4-currents J* is 

q{r) = r'^E{r) = An f J^r'^e^^+^'^'^dr, (5) 
Jo 

E being the intensity of the electric field. Here p and p are the pressure and matter- 
energy density, respectively. We have used prime to denote derivative with respect 
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to radial coordinate r o nly. The field equations without the charge q have been 
studied by Buchdahl and Bayinl^who found physically meaningful solutions. 
Equating equations ([2]) and ([3]) we get a differential equation by assuming = 
and B<^^ in the form 

$"-c— = - $''+1^, (6 

where $ = $(7") and c = (^fe^ - ^a^ - afe + 6 - a)/(6 - a). 

Hence, to solve the above equation ([6]) for $(r) we make use of the ansatz 
q{r)'^ = K^^{r)'^ so that the above equation becomes 

with c = b+l + d where d behaves as the polynomial index. A direct solution for the 
above second order differential equation is difficult to find out. However, we can set 
different conditions for the parameters to get a solvable equation. Let us therefore 
consider the following two cases. 



2.1. The case for c — 1 and d — —b 

For the specification c = 1, when d = —b, we get the following solutions 

2.2. The case forc^l and d = + a^)/{b - a) 

For this choice of c and d the solutions set becomes 

B{l-c) 



$1-= = 



A 



e'' = B 



2(6 -a) r2 
B{1 - c) 



2(6 -a) r2 
'B{1 - c) 



Cor' + Ci 



-a/(l-c) 



-,6/(l-c) 



+ Cor2 + Ci 



(8) 
(9) 

(10) 



(11) 



(12) 



(13) 



2(6 -a) r2 

We notice from the above two subcases that Co and Ci represent integration 
constants and their expressions will be determined in Section 4 in terms of m and 
R. 
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3. Boundary conditions 

The exterior field of a spherically symmetric static charged fiuid distribution de- 
scribed by the metric (1) in isotropic coordinates is the unique Reissner-Nordstrom 
solution 



2m q 



2m 



dr^ - r^{de^ + sin^9d(t>'^), (14) 



which by the radial coordinate transformation 



takes the form 



ds 



4^/2 4^/. 



,dt^ 



_2i 

4r'' 



m q 
2?^ 2? 



2r' 



1 



2r' 2r' 



Therefore, matching at r' = i? we get 

( 



MR) 



1 _ m I _2 



(15) 



{dr^ + r^dO^ + r^sin^ed(l)^). 



(16) 



(17) 



(18) 



where R is the boundary of the star. At the boundary pressure is zero. Using this 
we will get an expression for R in the next section. 



4. The study of the solutions 

4.1. The case for c — 1 and d — —b 

8np + 6bCo + f (s:^^ + 2Cor)2 + ^]e-^^)- (20) 

where :DW = 2(£l)^ + Cor2. 

In Fig. 1, Fig. 2 and Fig. 3 we plot the Snp, Snp and [q{r)]'^ against the parameter 
r for some fixed values of the integration constants Cq and Ci and parameters A, 
B, a, b and K. 

Fig. 4 above indicates that P + p is positive. That means though that the Weak 
Energy Condition (WEC) is violated in our case but obeys the Null Energy Con- 
dition (NEC) as well as the Strong Energy Condition (SEC) and the Dominant 
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C, = 1, C = 1, B = 0.5, A = 0.5, a: = 1, a = 1 



S Kp 



0.009- 

0.008 

0.007 

0.006 

0.005 

0.004 

0.003 

0.002 

0.001 

0.000 



b = 2 



Z> = 1.5 



Fig. 1. Bayin's isotropic solutions with charge showing pressure vs radius plot for the given 
specifications in the figure (Case-4.f). 



Energy Condition (DEC). This is a very specific charged fluid solution. In reference 
to the figure for p (Fig. 1), since p — gives the radius of the charged fluid (rather 
we would say, p{r = R) = 0), one may conclude that the radius of the charged 
fluid falls within 3 to 4 Km. Thus our charged fluid is highly compact and we may 
consider our charged fluid as typically a highly compact strange star (we call it 
'strange' in the sense that EOS components follow peculiar properties here i.e. they 
obey all energy conditions except WEC). In this context we would like to mention 
that in the original Bayin solutions, case V and case VI, p is also negative. 

Therefore, following Bayin El we are now in the position to find out pressure-to- 
density ratio at the centre of the charged sphere in the form 

El^J_^, (21) 
Pc b 

We notice that the pressure-to-density ratio at the centre of the charged sphere is 
entirely independent of the A and K parameters and integration constant Ci . In the 
case 6 — 2a > with 6 > we have ^ < and this is not the case with simultaneous 
positive pressure and density at the origin. For b — 2a < and 6 > we obtain 
^ > 0, which is a physically meaningful result for positive values for pressure and 
density at the origin. As a special case, if we choose a — 0, then we can easily recover 
the equation of state related to the Cases I and II for n = 1 of Ray et al. ^ which 
refers to the 'false vacuum' or 'degenerate vacuum' or 'p-vacuum' equation of state 
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Fig. 2. Bayin's isotropic solutions with charge showing density vs radius plot for the given spec- 
ifications in the figure (Case-4.1). 



Physical status a b b — 2a PcJ Pc 

Dust arbitrary 2a 

Vacuum fluid arbitrary b —1 

Stiff fluid 1 1-1+1 

Radiation 2 3 —11/3 



| i <^ | io | i!j | /U | .j-Q .j-jj^g repulsive nature of the pressure this provides a mechanism 
to avert the problem of singularity at the centre. On the other hand, if we choose 
a = 2 and 6 = 3 then we can recover the result of Ray et al. ^ related to the Case I 
for n = 3 which refers to radiation. In Table 1 we have shown different possibilities 
for different values of the parameters a and b. 

Substituting p{R) — 0, we can get an expression for the boundary of the charged 
sphere as follows: 

b{b - 2a)CiR^ + 2{b - a)Coi?« - '^'''ff f ^° i?^ + '^'^ff g'/' = 0. 

Though it is a complicated power law equation for i?, and hence very difficult to 
solve yet from this we get the same expression R^ — — ^^.^^^^a^Co ''^hi*^^ was obtained 
by Bayin (1978) in his non-charge case (see equation (4.23) there in). 

Again, comparing Reissner-Nordstrom metric in isotropic coordinates (16) [and 
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Fig. 3. Bayin's isotropic solutions showing charge vs radius plot for the given specifications in 
the figure (Case-4.1). 



hence (17) and (18)] we get 



(I + JIL)2 _ _2L 
\^ ^ 2Rj 



(22) 



so that 



where 



/ TO \ 2 




+ d - 


4i?2 



^ _ -G±(G2-4Fi/)i/2 
^0- 



F^b{b- 2a)R^ 



G = 2(6 - a)i?6 - ^JLJ^BK^R^ 

[b-a) 



(23) 

(24) 

(25) 
(26) 



(27) 
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Fig. 4. Bayin's isotropic solutions with charge showing EOS (pressure + density) vs radius plot 
for the given specifications in the figure (Case-4.1). 



and 



Ci = 



-2/a 



'1 _L ILL] _ r 



-2/a 



j^l/a^-D(R)^ 



Therefore, for B we get 



B 



1 + ^ 



1 



9^ ■ 


2(2o-6)/ 












^ -26/a 



(28) 



(29) 



Therefore, considering arbitrary values for a, h and A we get expressions for Co, 
Ci and B in terms of m and R. 

4.2. The case forc^l and d = -h(b^ + a^)/{b - a) 



b-2a /S'V b-a f S' 



4(1 -c)\S 



br \S 



(30) 



87rp = / 



5' 



4(1 -c)y V-S" 



2 

r I 5 



(31) 
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where 



B{l-c) ■ 



(32) 



(33) 



S = ; ^ + 26or, 



(34) 



5"' 



35(1 - c) 



2Co. 



(35) 



The pressure-to-density ratio at the centre of the charged sphere in this case is 
as fohows^ 



Pc 
Pc 



1. 



(36) 



T he same 'fa lse vacuum' or 'degenerate vacuum' or 'p-vacuum' equation of 
state | 17 | 18 | 19 |20| jg achieved again without imposing the condition a = in the 
present case. Therefore, also due to the repulsive nature of the pressure this pro- 
vides a mechanism to avert the problem of singularity at the centre. 

The radius of the star is given by the expression obtained from the condition 
p{R) — 0. However, here S and 5" have to be taken at r = i?, so that the final 
expression for the radius R can be written as 



(6 - 2a) 
4(1 -c) 



- I 



l-c 



-q 



(37) 



where we have used the symbol s = S{R). 

We can find out the values of Cq, Ci and B in terms of M and R from the 
equation (37) and the following boundary conditions 



1 



(38) 



Bs'/^ 



l-c) 



2r) 4i?2 



(39) 



Due to the complicated nature of the solutions above, we could not discuss these 
analytically, especially we are unable to find out the polytropic relat ion as obtained 
by Bayin ^ in his equation (4.43) for the classical limit Pc/ Pc 0^2II_ We expect 
to provide this study in a future project. 
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5. Conclusions 

The study of static charged fluid spheres in general relativity is far from being 
complete, but some recent studies have motivated us to extend one previous work'^ 
for improving the understanding of this topic. We found out the solutions of field 
equations in isotropic cases of charged fluid spheres in general relativity and per- 
formed a detailed study of the analytical solutions. In order to obtain the analytical 
solutions of the Einstein-Maxwell field equations, we made use of some assumptions 
and introduced two specialization, the first with the specification c = 1 and d = —b 
and the second with the values c ^ I and d — —^{b'^ + a'^)/{b — a). We calculated 
the corresponding expressions for the function (j){r) and for the metric coefficients 
e"^ and e^, which depend on the parameters A, B, a, b and integration constants Co 
and Ci- Matching the interior solution to the exterior Reissner- Nordstrom metric 
in isotropic coordinates at r = i? we obtained the formulae for e'^^^' and e^''^\ 

We performed a deeper investigation of the solutions and for the case c = 1 
and d = —6 we established the expression for the pressure-to-density ratio at the 
origin given by ^s. = —h^^liL and the imposed conditions for a and b for obtain- 
ing a physically meaningful model for the star. For a zero value of aweobtained 
the case of the false vacuum or degenerate vacuum or p-vacuum | 17 | 18 | 19 120] 
addition, we found out the explicit expression for the radius R and discussed the 
physically reasonable cases. The radius R for the charged sphere has a completed 
expression which, however, reduces to the same expression as in the non-charge case 
of Bayin Finally, comparing Reissner- Nordstrom metric in isotropic coordinates 
we established the expressions for the integration constants Cq and Ci and B pa- 
rameter in terms of mass m and the radius of the star R. On the other hand for 
the specification c 1 and d — — + a^)/(6 — a) we computed the pressure p 
and the density p that present a dependence on a new function S{r) and its first 
and second derivatives with respect to the r coordinate. In this case the expression 
for the pressure-to-density ratio at the origin can directly be given by y- = —1. We 
also gave the expressions that allow to evaluate the boundary R of the star and the 
values of Co, Ci and B in terms of m and R, but due to their complicated form we 
plan to perform a detailed study in a future work. 

The most straight forward observation in the present work is that from our 
results we are able to recover the neutral cases of Bayin We also observe that 
the models presented here have negative energy density and positive pressure and 
support the conclusion that static charged fluid spheres are connected to compact 
strange stars. We conclude that in the last years there is progress in the study 
of static charged fluid spheres which are connected to compact stars, but deeper 
and wider investigations for searching for new solutions need to be performed. One 
possibility will be to investigate our scenario based on different specialization and 
obtain new analytical solutions that can lead to other particular cases and more 
constrained connections with the involved parameters. 
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